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Abstract 

We compute the thermal free energy in large N U (N) Chern-Simons-matter theories 
with matter fields (scalars and/or fermions) in the fundamental representation, in the 
large temperature limit. We note that in these theories the eigenvalue distribution of the 
holonomy of the gauge field along the thermal circle does not localize even at very high 
temperatures, and this affects the computation significantly. We verify that our results 
are consistent with the conjectured dualities between Chern-Simons-matter theories with 
scalar fields and with fermion fields, as well as with the strong-weak coupling duality of 
the J\f = 2 supcrsymmctric Chern-Simons-matter theory. 
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1 Introduction 

An interesting class of conformal field theories in three dimensions comes from coupling a 
(topological) Chern-Simons (CS) theory to massless matter fields. Since the Chern-Simons 
coupling is quantized and cannot be renormalized beyond the one- loop level, it is easy to 
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make such theories conformal by tuning away any relevant operators. When the matter fields 
are purely fermionic this is enough to make the theory conformal. When scalar fields are 
present there are also classically marginal couplings (of the schematic form <^> 6 or (j) 2 iptp) and 
one has to arrange for their beta functions to vanish, but often this can be done, at least in 
the large N limit. 

In the past year, a specific class of conformal Chern-Simons-matter theories has been 
investigated in detail [1, 2, 3, 4, 5, 6]. This is the large N 't Hooft limit of a U(N) Chern- 
Simons theory at level k, coupled to a single matter field (scalar or fermion) in the funda- 
mental representation (the generalization to an O(N) theory and/or to Nf matter fields is 
straightforward). The large iV limit involves taking both N and k to infinity, keeping fixed 
the ratio A = N/k. These theories come in two versions; a "regular" version where one does 
not turn on any other couplings except the gauge coupling, and a "critical" version which 
can be defined by turning on a quartic coupling and sending its coefficient to infinity, or by 
adding a mass parameter and making it into a dynamical field that is integrated over. It was 
found that this class of theories has some very nice properties : 

• For A = the theory is free and has an exact high-spin symmetry, with spin s generators 
for every positive integer s. For finite A and large ./V this symmetry is still approximately 
present, in the sense that the anomalous dimensions of the high-spin currents are of order 
1/JV [1, 2]. 

• The approximate high-spin symmetry allows for an exact evaluation of all 2-point and 
3-point functions of "single-trace" primary operators 1 , up to two (or three) parameters 
[7, 8]. 

• The 2-point and 3-point correlators of these theories can also be computed exactly in 
the planar limit by resumming Feynman diagrams, and agree with the general results 
based on symmetries [5, 6]. 

• It was proposed [1, 2] that these theories (as well as various supersymmetric extensions 

[3]) have a gravity dual given by a parity breaking version of Vasiliev's higher-spin 

theory on AdS± [9, 10], with a parity-breaking parameter 6 that is a function of the 

't Hooft coupling A (the exact gauge theory computation of [5], direct 3-point function 

calculations in the bulk [11] and arguments in [3] suggest a simple linear relation between 

the bulk phase 6 and A). This duality is a generalization of the original conjectures 

[12, 13] relating the singlet sector of bosonic/fermionic vector models to the parity 

invariant Vasiliev's theories. 

1 We will refer to gauge-invariant operators coming from the contraction of a single fundamental and a 
single anti-fundamental field (with any number of covariant derivatives) as single-trace operators; their large 
N properties are analogous to those of single-trace operators in theories of adjoint fields. 
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• The gravity dual is the same for the theory of scalars coupled to Chern-Simons, and 
for fermions coupled to Chern-Simons, suggesting a duality between these theories, and 
this duality is consistent also with their planar 2-point and 3-point functions. More 
precisely, the "critical" version of the scalar theory has the same gravity dual as the 
"regular" version of the fermionic theory, and vice versa. 

In the large ./V limit, the considerations above suggest an exact duality between the 
"critical" U(N)k YM CS theory coupled to a scalar field, and the "regular" U(kyu)-N CS 
theory coupled to a fermionic field; here we are using the definition of the coupling constant 
kyM that arises when the theory is regularized using a Yang-Mills theory at high energies 
(this is the definition that agrees with the level of the affine Lie algebra that arises on 
boundaries of the Chern-Simons theory, but many formulas depend on the renormalized 
coupling k = kyM + N sign(A;YM); in this paper we will generally use this renormalized 
coupling in all formulas, unless indicated otherwise). In the absence of the matter fields this 
duality is simply the level-rank duality of pure Chern-Simons theories [14, 15, 16, 17, 18, 19] 
(and, whenever they are non-zero, the Chern-Simons contributions dominate in the large iV 
limit over those of the matter fields), but the claim is that there is still an exact duality also 
after adding the matter fields. The duality was conjectured in [5] to hold also for finite values 
of N, with a small shift in the CS level on the fermionic side. 

The tests of this (non-supersymmetric) duality so far involve exact computations of var- 
ious 2-point and 3-point correlation functions in the planar limit, which agree precisely. 
However, as shown in [8] , these correlators are all determined by the almost-conserved high- 
spin symmetry, and it would be nice to be able to test the duality by properties that are not 
determined purely by this symmetry. This is particularly important since the gravity dual 
Vasiliev theory appears to contain additional coupling constants that could affect 5-point and 
higher correlators of single-trace operators, and their role (if any) in these dualities is still not 
clear, so one may wonder whether the two theories may have the same 2-point and 3-point 
functions for their single-trace operators (determined by the high-spin symmetry) but still 
not be identical. 

In this paper we provide further evidence for this duality by matching the thermal free 
energy between the two theories mentioned above, at large enough temperatures and in the 
planar limit. The free energy for the fermionic theory was computed already in [ ], while the 
scalar computation (by the same methods) was performed in [1, 20]. These computations do 
not agree with the duality (nor with supersymmetric versions of the duality, for which there 
is much more evidence); the scalar computation does not even show a vanishing free energy 
in the strong coupling limit A — > 1 where (T(x)T(y)) goes to zero, indicating that it cannot 
possibly be correct. 

The computations in the literature assumed that the eigenvalue distribution of the holon- 
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omy of the gauge field around the thermal circle localizes around zero at high enough tem- 
peratures, so that the periodicity of the fermionic/scalar fields can be taken to be the naive 
one (without any shift coming from the gauge holonomy). For free theories this is certainly 
correct [21, 22]; the matter fields tend to cause the eigenvalue distribution to localize, and 
while on some spatial manifolds (like S 2 ) there is a repulsive term coming from the measure, 
the effects of this term go to zero at high temperatures. For asymptotically free gauge theories 
it is also clearly true that the holonomy eigenvalues localize at high enough temperatures. 
However, for conformal theories at finite coupling it is not clear if this is correct or not. 

In section 2 we argue that the correct large N eigenvalue distribution for the thermal gauge 
holonomy in a Chern-Simons-matter theory of the type described above is not localized, but 
actually is a uniform distribution with a width covering a fraction |A| of the unit circle. 
Recall that in our convention A goes from zero at weak coupling to one at strong coupling, 
so that at weak coupling the eigenvalues do localize as expected, while in the strong coupling 
limit the eigenvalue distribution becomes uniform (as in a confined phase, though there is 
no confinement here). This effect comes from the Chern-Simons interactions, that make the 
eigenvalues effectively fermionic so that they cannot be too dense, despite the attractive forces 
coming from the matter sector. 

In sections 3-5 we compute the thermal free energy for the scalar and fermion theories, 
both "regular" and "critical", using this eigenvalue distribution. The computation is com- 
pletely analogous to the computations of [1, 4, 20], just with a different periodicity condition 
for the matter fields, arising from the holonomy. We show that with this correct eigenvalue 
distribution the free energies have the expected properties (going to zero in the strong cou- 
pling limit), and that there is a precise matching between the thermal free energies of the 
scalar and fermion theories. This provides an additional test for the duality between these 
two theories, that does not follow just from their approximate high-spin symmetry (as far as 
we know). 

In section 6, we then generalize the computation to a theory that has both scalars and 
fermions, with general couplings. In particular, we verify that in the special case of the theory 
with J\f = 2 supersymmetry, the thermal free energy in the large N limit is invariant under 
the supersymmetric version of the duality [23, 24]. 

Since this paper is rather long, we summarize our results for the free energies of the various 
theories in section 7. Our analysis of the critical theories automatically gives us the thermal 
free energies also for the Chern-Simons theories coupled to massive scalars and fermions, and 
in this section we also briefly discuss some of their properties. 

Our methods for computing the holonomy apply to any Chern-Simons-matter theory in 
which the number of matter fields is much smaller than N 2 in the large N limit. It would be 
interesting to understand their implications for theories like those of [25] where the number of 
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matter fields is of order N 2 , so that one cannot neglect their back- reaction on the eigenvalue 
distribution. This should be important in any computation of the (planar) thermal free 
energy in these theories. Our methods are sufficient for computing the order A 2 correction to 
the free energy in these theories. 

It would be interesting to understand how to reproduce our results for the thermal free 
energy (or various other effects of the Chern-Simons part of the conformal field theory [26, 27]) 
in the gravity dual theory. It would also be interesting to find additional tests for the field 
theory duality, in particular for finite values of N. 

2 The Holonomy Along the Thermal Circle 

When one studies the thermodynamics of gauge theories it is common to have long distance 
modes that require special treatment [28, 29, 30]. These are captured cleanly by thinking 
of the thermal computation as a Kaluza-Klein reduction of the Euclidean theory on a circle. 
If we go to energies lower than the inverse radius of the circle (the temperature), we might 
encounter light degrees of freedom that contribute to the free energy. In particular, the 
holonomy of the gauge field on the circle is classically a massless field. These light degrees of 
freedom need to be taken into account correctly. 

In our case, we are interested in conformal field theories coming from large N U(N) 
Chern-Simons theories coupled to matter fields in the fundamental representation. On general 
grounds, the free energy of a 3d conformal field theory at temperature T and volume V is 
proportional, for large enough volumes, to VT 2 . For a standard gauge theory coupled to 
fields in the fundamental representation, the large N free energy would be dominated by the 
contribution of the gauge fields (scaling as -/V 2 ), while the matter field contributions scale as 
N. In our case, the Chern-Simons theory is purely topological, so the contributions of this 
theory on its own are independent of the volume, and one thus expects the free energy to 
scale as N. One may expect that at large enough volumes (or, equivalently, large enough 
temperatures) the contribution of the matter fields will be dominant, and will make the 
eigenvalues of the holonomy matrix of the gauge field localize near zero (as in asymptotically 
free gauge theories). This is the assumption that was made in previous computations of the 
free energy. However, we will see that the presence of the Chern-Simons theory modifies this, 
even at very large volumes. 

In order to analyze the thermal partition function, we reduce the Euclidean theory on 
the thermal circle of circumference f3 = 1/T. For the purposes of this section, we begin by 
integrating out the matter fields in the fundamental representation (at finite temperature 
and coupling we expect all these fields to be massive) . We are then left with an effective two 
dimensional theory, involving the dimensional reduction of the gauge field. We have the two 



5 



dimensional gauge fields, and also the holonomy of the gauge field along the thermal circle, 

dx 3 A 3 . (2.1) 

The i is introduced so that a is hermitian (we use a convention in which the fields are 
anti-Hermitian) . 

Let us first consider this dimensional reduction for the pure Euclidean Chern-Simons 
theory, 

S ° S = trJ Tr ( A ^ A P ~ J A » A " A p) . k e 2 • (2-2) 
In that case we get an effective action of the form 

S = — [ d 2 xTr(Fa), (2.3) 
2vr J 

where a is the holonomy (which is a scalar field in two dimensions), and F is the field strength 
along the two non-compact dimensions. Since the original action is topological, this action is 
also topological. 

Let us now compactify one of the IR 2 directions on a circle of circumference L, and think 
of the other non-compact direction as (Euclidean) time. We now have a second holonomy, 

: f dx 1 A 1 . (2.4) 
J o 



b = i 



We can gauge fix the other component of the gauge field Aq to zero, and the two dimensional 
action then becomes a quantum mechanics action for the two holonomies, of the form 

S = — fdtTr(ad b). (2.5) 
2vr J 

This has the form of a matrix quantum mechanics action with a single matrix degree of 
freedom, since b is canonically conjugate to a. We can diagonalize a and consider only the 
dynamics of the eigenvalues. As is standard in matrix quantum mechanics, this introduces a 
Vandermonde determinant from the measure, which we can swallow into the wavefunction, 
but this makes the remaining part of the wavefunction anti-symmetric in the eigenvalues, so 
they need to be treated as identical (due to the remaining Weyl symmetry) fermions [31]. 
Thus, we get a Euclidean action for N fermionic variables [32], 

k N r 

s=yH dthiQi - (2 ' 6) 
n i=i 

For each eigenvalue we have the action for a particle in a magnetic field, on the lowest 
Landau level. Due to the fact that a, and bi come from holonomies, they all have periodic 



6 



identifications, which in our normalization are a% ~ a% + 2ir and 6, ~ 6j + 2tt. So we have 
iV fermions moving on a torus, with a magnetic field that has k units of flux. Having flux 
k implies that each Landau level has precisely |fc| states. So, we have N fermions filling 
states. In our conventions k is the renormalized Chern-Simons coupling (which for the 
SU(N) part of the gauge group is related to the Chern-Simons coupling that arises from a 
Yang-Mills regularization by k = fcyM+^sign^YM)), so we always have \k\ > N, and we have 
(!*!) different ways of filling these states 2 . In the pure Chern-Simons theory the Hamiltonian 
vanishes, so all these states have the same energy (zero energy), giving the ground state 
degeneracy of the U(N) Chern-Simons theory on a spatial torus [32] . 

Up to now our analysis was exact. Next, we add the matter in the fundamental repre- 
sentation of the gauge group, in the large N limit. We imagine that the thermal circle is 
much smaller than the other circle, and we expect all the matter fields to obtain masses of 
the order of the temperature. When we integrate out the matter fields we get a contribution 
to the free energy of order N times the volume. In this computation (which we will perform 
in detail in the subsequent sections for a specific holonomy) we treat the holonomy as a con- 
stant background field, so we obtain a result which depends on the holonomy. In the limit 
where the thermal circle is much smaller than the other circles, the other effects of the gauge 
field are very small. The most important term that we get in the effective action is then a 
potential for the holonomy 



where V{a) is some (gauge-invariant) functional of the holonomy, which scales as N in the 
large N limit. (Additional terms depending on the two dimensional gauge fields also arise 
when integrating out the matter fields, and we will discuss them below.) At least to leading 
order in the 1/k expansion, the potential V(a) has a minimum at a = 0, and the expansion 
around zero looks like V(a) — V(0) oc Tr(a 2 )T 2 . 

One may be tempted to set a = 0, as was done in previous computations in thermal Chern- 
Simons-matter theories. However, recall that when we consider this theory compactified on 
an extra S , we saw that the eigenvalues of the holonomy a behave like the positions of 
fermions in the lowest Landau level. Thus, we cannot set them all to zero; the state a% = 
is not consistent with the Pauli exclusion principle. The best we can do is to fill a Fermi sea 
so that the N levels are as close to zero (in the a direction) as possible. Since kai/2ir is the 
conjugate momentum to 6j, which has periodicity 2ir, if we work in a basis of eigenstates of 
aj, then the eigenvalues are quantized in units of (27r/|fc|). We want to fill the N eigenstates 
that are closest to a = 0. In the "phase space" torus parameterized by a and b, this means 

2 Note that there are other regularizations in which one gets a similar action (2.6), but with a coefficient kyM 
(for the SU(N) part) and with bosonic eigenvalues [33]. Note that there are no Vandermonde determinant 
factors for the Chern-Simons theory on T 3 [34], but the eigenvalues still behave as fermions in this case [32]. 
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that we are filling all states in a band of width 2-7r|A| around a = 0, and uniform in the b 
direction, see figure 1(b). Notice that in the large N limit, we have many eigenvalues and we 
can view the distribution of fermions as a fluid, ignoring the details of the wavefunctions of 
these fermionic eigenvalues on the torus. 




a 



b b 

(a) (b) 

Figure 1 : (a) One out of many possible configurations of fermions in the lowest Landau level 
on the torus, (b) The potential generated by integrating out the matter fields selects the 
configuration where the fermions are localized in a band around a = 0. The fact that L S> /3 
breaks the symmetry between the two holonomies. 

The matter contributions break the degeneracy between the different ground states of the 
Chern-Simons theory, and imply that this specific state has the lowest energy. We will argue 
below that the other contributions from integrating out the matter fields do not affect this 
computation at leading order in the large N limit. In principle we should now sum over all 
the different states with a thermal distribution. The energy density difference between the 
lowest energy state that we found and the next state is of order T 3 /\k\ (since we can shift a 
single eigenvalue by 2ir/k). Thus, if VT 2 is much larger than N and k, we can focus on the 
ground state, and ignore all other contributions to the free energy. From here on we will be 
working in this limit for simplicity. The regime VT 2 ~ 1 was discussed in [26, 27]. 

The discussion above implies that at large volumes/temperatures we should diagonalize 
the thermal holonomy and distribute the eigenvalues uniformly in a band of width 27r|A| 
around zero. It is important that even though to derive this result we compactified an extra 
circle, this width is independent of the radius of the second spatial circle that we used for 
this derivation, as long as L > So, this distribution of holonomy eigenvalues affects the 
extensive part of the free energy (the part that is proportional to the radius of the second 
circle). It turns out that this change in the eigenvalue distribution affects the thermal free 
energy already at order A 2 in the weak coupling expansion. As a consistency check on this 
value for the holonomy, one can check that with this eigenvalue distribution the traces of the 
holonomy matrix e a transform in the correct way under the level-rank duality (the traces 
of this holonomy in symmetric representations in the theory with coupling A are mapped 
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to the traces of the holonomy in anti-symmetric representations in the theory with coupling 
A — sign(A)). This would not be the case for a = 0. 

It seems surprising that a finite-volume effect like the eigenvalues being fermions modifies 
the free energy even in the large volume limit. We can obtain an alternative perspective on 
this by following the standard treatment of long distance modes in thermal gauge theories 
[30], adapted for Chern-Simons matter theories. Again, we consider for simplicity the large k 
limit (with large N, and with A = N/k fixed and small). Then, after dimensionally reducing 
on the circle we get the following effective action 



S = I d 2 x 



^Tr(Fa) + T 2 Tr(a 2 ) + 

Z7T 



(2. 



where we expanded around a = and kept just the leading terms, ignoring coefficients of 
order one. We can now integrate out the scalar field a, obtaining 



<7 k2 

rp2 



jd 2 xTi{F 2 ). (2.9) 



This is a two dimensional Yang-Mills action with an effective coupling g\ D oc So for 
large k this is very weakly coupled at the scale T. In two dimensional Yang-Mills there are 
no propagating degrees of freedom, but there are non-trivial flux excitations. Here we are 
interested in the ground state energy of this effective two dimensional Yang-Mills theory. We 
again compactify the x 1 direction on a circle of circumference L. Choosing the other direction 
to be time, we end up with a matrix quantum mechanics for the holonomy b = i f dx 1 Ai 
around the x 1 direction. Again we have free fermions by diagonalizing the matrix; this time 
we have free fermions on a circle. The effective Lagrangian takes the form 

r i P N 

Socl dt-^Y< b l ( 2 - 10 ) 

i=i 

Since b{ is periodic with period 27r, the quantized momenta for 6, have the form = with 
integ er Tij, where pi ~ b% . We fill the Fermi sea up to \n\ — N/2, and we end up with a 
ground state energy of the form 

^ K E ^r n f^^^^ LT2Nx2 - ( 2 - n ) 

m=-N/2 

In particular, we see that it is proportional to L, so it is an extensive contribution. And, we 
see that it is (at leading order in perturbation theory) an order A 2 correction to the leading 
order in N result for the free energy (that would come from the vacuum energy of the matter 
fields themselves); thus, it cannot be ignored in the large limit. As usual, in order to speak 
in a meaningful way about the ground state energy we need to define a subtraction scheme. 
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The fact that the scheme we used here is the correct one is most transparent according to our 
previous discussion, where we started from the dimensional reduction of Chern-Simons theory 
and viewed the same correction to the energy as arising from integrating out the matter fields 
in the fundamental representation. 

As we integrate out matter fields, we also induce terms of the form 



with order one coefficients (depending on A) and with powers of the temperature fixed by 
dimensional analysis. Such terms do not affect the previous discussions, because they are 
subleading in the 1/N expansion. 

The spreading of the eigenvalues of the holonomy is expected to be a general feature in the 
thermodynamics of Chern-Simons-matter theories. For a finite number of matter fields in the 
fundamental representation we saw that the eigenvalues have a uniform distribution with a 
width proportional to |A|, with corrections that are small at large N. It is more complicated 
to analyze what happens when the number of matter fields is of order N 2 ; in that case 
corrections like (2.12) are no longer negligible, and could affect the eigenvalue distribution. 
The precise eigenvalue distribution could also depend on the precise field content, since (for 
instance) for adjoint matter fields the potential V(a) depends only on differences between 
eigenvalues. In any case, we expect that there should be a non-trivial eigenvalue distibution 
also in such cases, which would give extra terms starting at order A 2 in the thermal partition 
function. In particular, these effects should be added to the results in [35], and to any other 
perturbative computations of the thermodynamics of Chern-Simons-matter theories. 

The derivation we have presented here is most clear if we consider the free energy for the 
theory on T 2 x So. In this case the T 2 gives us the extra circle that we needed to run our 
argument, and we can take the large volume limit of the T 2 and obtain the extensive part of 
the free energy as we discussed above. One could wonder if similar effects also appear when we 
consider the three dimensional theory on £2 x S\ for other Riemann surfaces £2, m the limit 
that Sp is very small. Pure three dimensional Chern-Simons theory on such manifolds was 
studied in detail in [34]. They argued that the theory can be reduced to a two dimensional 
theory of the form (2.3), where now the integral is over £2- In. addition, we can further 
diagonalize the holonomy in the circle direction. The resulting effective action is independent 
of the manifold, up to overall factors (involving powers of the Vandermonde determinant) 
which do not scale like the volume of £2, and thus do not contribute to the free energy in 
the limit of large £2 volume. Since the analysis on the torus gave us an extensive effect, we 
expect that the same effect should arise on a general Riemann surface £2- It would be nice to 
show in detail how this works. A particularly interesting case is when £2 = S 2 . It was argued 
in [36] (for A = 0) that the thermal circle needs to be very small, of order 1/y/N as compared 
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to the size of the S 2 , in order to obtain an answer for the free energy that is linear in N. Note 
that here we are also discussing very large temperatures of the same order. On the sphere, 
there is a volume independent contribution to the effective action for the thermal holonomy, 
a, that leads to a repulsion between eigenvalues. This will lead to an eigenvalue distribution 
that is uniformly spread over the whole circle for small temperatures, but for very large 
temperatures this is overwhelmed by the effects we encountered on the torus. For negatively 
curved S2 the extra volume independent factors lead to attractive interactions between the 
eigenvalues, so we expect the eigenvalues to localize already at smaller temperatures. 

In the remainder of this paper we will compute the resulting free energy exactly as a 
function of A, but to leading order in the 1/N expansion. We work on R 2 , which can be 
thought of as the large volume limit of T 2 (as we mentioned, we expect the result to be the 
same for large enough volumes on any two dimensional space). The computation is similar 
to those done in [1, 4, 20]. We simply need to repeat all their steps, but now with a constant 
holonomy a matrix, with the eigenvalue distribution that we found, in which the eigenvalues 
of the holonomy are distributed in a band as in figure 1 (b). We view this as a constant 
background field in the perturbative computation of the free energy. 

3 Preliminaries 

In the following sections we consider the theory of N complex scalars or Dirac fermions, in the 
fundamental representation of U(N), coupled to a U(N) gauge field Ap with Chern-Simons 
interactions in 3 Euclidean dimensions. The Chern-Simons action (2.2) can be written as 3 

Scs = ~j d 3 x (Ald v A a p + l -f abc AlA b v A c )j , fceZ. (3.1) 

We take the 't Hooft large ./V limit, keeping A = £ fixed, and work (following [ ]) in light-cone 
gauge in the x 1 — x 2 directions, A- = 0. 4 In this gauge the gluon self-interaction vanishes, 
and the Chern-Simons action reduces to 

Scs = ^ Jd 3 xA a + d.A a 3 . (3.2) 

We will regularize our theory by placing a sharp momentum cutoff A in the x 1 — x 2 directions. 
This will cure most UV divergences, and the remaining ones (that will occur when computing 
determinants) will be handled as they are encountered. In this regularization (used also in 
[1, 2, 4, 5]), k is related to the level &ym that appears in a regularization by a Yang-Mills 

3 The generators T a are anti-Hermitian, with the normalization Tv(T a T b ) = ~\5 ah . 

4 Light-cone coordinates are defined by x = (x 1 ± ix 2 )/\/2. The invariant tensors are given by #33 = 
5+_ = 1 and e_| 3 = i. 
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action (and that appears, for instance, in the Wess-Zumino-Witten model on boundaries) by 
k = /cym + ^ r sign(/cYM) (for the SU(N) part). 

We will take our manifold to be R 2 x S , where the compactified direction is x 3 and 
has length j3. As argued in Section 2, we should take the gauge field to have non-trivial 
holonomy around S* 1 (the eigenvalues of the holonomy are gauge-invariant). We can describe 
the holonomy in terms of the zero mode 

A Z = TF I d 2 xA 3 , (3.3) 



V 2 



where V2 is the regularized volume of R 2 . Using a residual gauge transformation (that pre- 
serves the light-cone gauge) we set 83 A3 = and diagonalize .A3. In the notations of Section 
2, the holonomy along the thermal circle is 

a = i(3A 3 , (3.4) 

whose eigenvalues are real and live on a circle of length 2ir, due to the periodic identifications 
which follow from the gauge symmetry. 

Most of the calculation can be done for general fixed a, but in the end we are interested 
in the large iV limit where the eigenvalues of a are uniformly spread around with width 
27r I A I , namely we take 



a(u) = 27r|A|u , u G 



1 1 

2' 2 



(3.5) 



and replace sums over eigenvalues by integrals over u with uniform distribution over the 
spread of eigenvalues. In other words, we replace 



N „J 



i=l " 2 



2 



f(a u ) — > N / dn/(2vr|A|n) . (3.6) 



1 



4 Bosonic Theories 

In this section we consider (following [2, 5]) the theory of N complex scalars <f>, coupled to 
the gauge field A^, defined by the action 

S = S sca i ar + S cs , (4.1) 

Scalar = J d 3 x (D^(D^) + ^(^4>f , = d„ + . (4.2) 

We include the coupling A6 which is exactly marginal in the large N limit [2]. We treat the 
holonomy as part of the free theory. Its effect is to shift the scalar momentum. Let us define 

P l i=P/i- i-Af, , (4.3) 
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where ^.3 is the zero-mode (3.3) and A± = 0. We will diagonalize the matrix ^.3 so that we 
have p3 = P3 — an/ /3 for the z'th scalar. The Feynman rules for this theory in momentum 
space include (in the A_ = gauge) 



j =(r 2 h 




= G Ufl (p)5 ab 



= i(T%+p f ,T a ) tJ 



= {T a ,T b }ijS tl 36 V 3 



i 3 



The non-vanishing components of the gluon propagator are 

4-7T? 1 

G +3 (p) = -G 3+ (p) = — — . (4.4) 

K p + 

It is important that this gluon propagator does not depend on p^. Thus, it is not affected by 
the holonomy. 



4.1 Exact propagator 

The exact scalar propagator is 5 

(0(p)fyt(-g)i> = 

where (S^)jj is the sum of 1PI diagrams. As we will see, at large ./V all contributions to 
(Ss)ji are proportional to the unit matrix in color space, so £#(p)jj = T,B(p)5ji. Also, 
none of the contributions will depend on p$. Using this and the fact that X# must be 
invariant under rotations in the x 1 — x 2 plane, we see that £5 can depend on p only through 
Ps = (Pi) 2 + (P2) 2 = 2p + p~ ■ Using dimensional analysis, we can write 

Sb(p) = -/3" 2 /u|(A, A 6 , p Ps , a) , (4.6) 
5 For the thermal theory P3 = 27rn//3, q = lim' //3, and we should interpret 2%S(p3 — 53) — > fi5 n:n i . 



P 1 ~ 



(2ir) 6 5 6 (p-q), 



(4.5) 



1 1. 
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and hbP 1 = ^bT is the thermal mass of the scalar. 6 

The bootstrap equation for is shown in figure 2. The 1-loop diagram (not shown since 
it vanishes) can be written as 

1 

' q 2 - £ B (<?) 

It is understood here and below that the full integration measure on R 2 x S 1 is 

1 f d 2 q 



V 3 q 



{T a q u + p u T a )- 2 ^^{T a p fl + tf T a ) 



(4.7) 



V 6 q 



(4.t 



where momenta q on S 1 are quantized as q n = ^p. It is easy to check that the diagram (4.7) 
vanishes using the relation for generators in the fundamental representation of U(N), 

1 



ij 1 kl 



and the fact that £ B and p are diagonal in color space for any vector p. 



■0- 




+ 



(4.9) 



(c) 



+ 



(d) 



Figure 2: Bootstrap equation for the scalar self-energy. A filled circle denotes the full scalar 
propagator, (d) corresponds to the insertion of the Ag(^ 2 ) 3 interaction vertex. 

In the remaining diagrams, the (f)du(jr vertex always appears with [i ^ 3, and therefore 
the Feynman rule for the vertex becomes simply iT^(p + p')^. The remaining diagrams can 
be written as 



(a) 



k J J q (l-p) + {q~p) + 



rpb 



I {T',T b }- 



w = («) = (?) 



rpa _ 



l 2 -z B (l) 
1 



q 2 -Z B (q) 



1 



(d) 



l 2 -Z B (l) q 2 -Z B (q) 
1 



{T a ,T b } 



j i 



J2«2 



q 2 -^ B (q)\ nii Yl 2 -T, B (l) 
[&ij&i-Lji&nj2 + (5 permutations)] . 
6 From now on, we may also refer to the dimensionless quantity /j,b as the "thermal mass". 



(4.10) 
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In diagram (d) the one-half is a symmetry factor. For general matrices A, B we have, using 
(4.9), 



T b A{T a ,T b }BT a = T b AT a B{T a ,T b } = -[BA + Tr(A)Tr(B)l] . 



(4.11) 



Using this in diagrams (a)-(c), we see that the first term of (4.11) is subleading at large N 
and can be dropped. In diagram (d), only the first permutation survives at large N. The 
diagrams can therefore be written as 



(a) 



47T% 



1 



(4.12) 



(b) = (c) = <\„ / V'HV'q 



/' 



3^3 (i + p) + {q + iy 



4V 

.] " ' 1 (l-p)+{q-l)+ 



Tr 



Tr 



P-E B (l)J ' 



(d) 



1 

iV 



D^Tr 



1 



We can simplify things by taking the derivative of the bootstrap equation 

£b(p)<% = (a) + (b) + (c) + (d) 
with respect to p~ , using the relations 

dp s p + 



d 1 



2^<5 z (p) 



dp p s 



(4.13) 



(4.14) 



The resulting equation is <9 p - £#(£>) = 0, and therefore and p? B do not depend on p s . 
Using this observation, we can simply evaluate the various contributions to Eg by setting 
p s = 0. With this choice, one finds that (b)=(c)=0 for symmetry reasons, and 



(a) = -A7r 2 X 2 5 ji 
A 6r 



1/2 
-1/2 



1 



(d) 



1/2 
1/2 



g 2 + /3- 2 /i| 

2 



du / £> 3 g- 



1 



(4.15) 
(4.16) 



where we have replaced the trace with an integral over the uniform spread of eigenvalues. 
The sum over the q% momentum takes the form 



^ oo 

F (<ls,U) = - 727 



n=— oo \ /3 



(f (n +|A|«))2 + g2 +j8 -2 M 2 



Re coth 



2j(/3g s ) 2 + /4 



(/3(? s ) 2 + //| + 27ri|A| 



(4.17) 
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The bootstrap equation for the self-energy can then be written as 

2 



S B = —A 



du I dq s q s F(q s ,u) 



where we defined 



1/2 
-1/2 



8"7T Z 



The thermal mass \xb then obeys the bootstrap equation 



|/xb(A, A 6 )| 



^ /-1/2 /-oo 
1/2 Jo 



Re coth 



X 2 + //^ 



x 2 + fj? B + 27ri|A|it 



(4.18) 



(4.19) 



(4.20) 



In writing this we made a change of variables x = f3q s , and took the square root of the 
equation for Eg- Using a cutoff A' = /3A, we find that the radial integral gives (using 
Re(log(z)) = log|z|) 



2 log 



sinh 



'x 2 + n 2 B + 27ri|A|n 



A' 



i=0 



A' - 2 log 



2 sinh 



\[j,b\ + 27ri| X\u 



+ o 



1 

(4.21) 



Subtracting the divergence using a mass counter-term, and taking the limit A — > oo, the 
bootstrap equation for /j,b now reads 



I MB | 







A / dulog 


2 sinh ^ 


J-l/2 





l/i^l + 27ri|A|u 



(4.22) 



This integral can be computed analytically in terms of dilogarithms, and we find 7 the 



equation 



±MB 



|Amb| 1 |A| 
2 2vri |A| 



Li 2 ( e -l^hHA|) _ L i 2 ( e - 



-|Msl+T«|Ah 



(4.23) 



In some cases there may be more than one solution to this equation (with different choices 
of the relative sign), and in such cases we should choose the solution that minimizes the free 
energy. We can always take (without loss of generality) ^b > 0, and we will do this in the 
rest of this section. The solution that minimizes the free energy is then always the one with 
the plus sign in (4.23). 



7 This is easily done by again replacing log \z\ — Re(log(z)). 

8 In this paper the dilogarithm function is defined on the principal branch, with a cut from 1 to +oo along 
the real axis. 
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To summarize, the exact scalar propagator is given by 



P 2 + /3- 2 /4(A,A 6 ) 



(27rf5 3 (p-q), 



(4.24) 



where the thermal mass fis is found by solving (4.23). Note that if we drop the integral over 
the holonomy in (4.22) and set u = 0, the equation determining the thermal mass agrees with 
the result derived in [4, 20]. 



4.2 Free energy 



It was shown in [4, 20] that the free energy has a simple expression in terms of the self-energy. 
In our case the calculation follows through, except that we should include the integral over 
the holonomy. At large N, we have 



^1/2 r j2 00 



log if - E B ) + 



2 E 



B 



3 p 



»2 



(4.25) 



2irn 



and p l j,=p ll w~h,n- 



where p 3 - — a ^ ^ M - Ffl j- 

Consider first the log term. The sum over n diverges, and we regulate it by introducing 
a negative contribution with a large mass /3 _1 M, so that we need to compute (here we use 
that E B = -/3- 2 n 2 B ) 



d 2 p 



p 2 + 

p 2 + p~ 2 M 2 



d 2 p 



io § n 



i+ 



Pi + W B 



2tt 



n — \\\u] 



(/i B -> M) 



d 2 p 

W) 



1 



2 log 



A' 



cosh(/3Jp2 + ^-2^) _ CO s(27t|A|«) 



(fJ-B -> M) 



2cosh(A/x 2 + /x|) - 2cos(27r|A|u) 



Oifl -»• M) 



(4.26) 



where again x = (3p s and A' = /3A. The first term has a UV power divergence, which we now 
subtract. We can write this term as 



1 



A' 



dxxl J x 2 + /i| + log \l + e-V^+^B - 2e-V^ 2 +M| C os(27r|A|u) 



(4.27) 



The divergence is now contained in the first term. The last term, which is finite, can be 



simplified by taking y = \/x 2 + ^l 2 b and sending A — > oo. After subtracting the divergence 
we find 



1 

27T/3 2 



^+ I dy y log (1 + e~ 2y - 2e~ y cos{2tt\ X\u)) 



(4.28) 
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1 /2 

Let us now compute the integral ; 2 du of this expression. With some work, one finds the 



result 9 



27T/3 2 



7Ti\X\ 



dyy 



l-'B 



Li 2 (e- y+7ri|A| ) - Li 2 (e" 



(4.30) 



This is the first term in (4.26). The second, regulator term in (4.26) gives a pure UV diver- 
gence (when taking M, A — > oo) that we subtract. This removes the determinant regulator. 

Let us now consider the second term in (4.25). Here the computation is similar to what 
we encountered in section 4.1. The bootstrap equation for the self-energy implies that 



Hb = ±27r|A| 



1/2 
-1/2 



du 



d 2 p 



E 



(4.31) 



and so we find that 



2 f 1 / 2 

-1/2 



du 



E 



d 2 p 



(2vr)2 ^ p2_ s 

v ' n=— oo 



Mb 



3vr/? 2 |A| 



(4.32) 



Combining (4.30) and (4.32), we find that the free energy is given by 



NV 2 j/4 
2vr/3 2 I 3 



1± 



+ 



\X\J ™\\\ 



dyy 



l-'B 



Li 2 ( e ^+HA|) -u 2 {e~ y - n i|A| )l \ , (4.33) 



where the choice of sign is correlated to the choice of sign in (4.23). Note that (4.23) is always 
given by setting to zero the /is-derivative of (4.33). Using (4.23) we can rewrite (4.33) as 



(3F B 



N V-2 



A 2iT 2 i(3 2 \ 3 
+ 



Li 2 {e~^ B+niX ) - Li 2 (e-^-" A ; 



dyy 



l-'B 



Li 2 ( e -w+^) - Li 2 (e- 



(4.34) 



It is obvious from this expression that the free energy vanishes (for fixed N) in the strong 
coupling limit |A| — > 1. The free energy with A$ = is shown in figure 4 below. 



4.3 The critical theory 

In this subsection we consider the critical bosonic theory, which at large N can be defined 
as the Legendre transform of the regular theory (with \q = 0) with respect to the scalar 
operator (jr <j>. In order to carry out the transform, we begin by computing the free energy in 

9 The following relation is useful here: 

Li 2 (e 2 ™) + Li 2 (V 2 "") = 2ti 2 B 2 {x) , < Re (x) < 1 , (4.29) 
where B2 (x) — x 2 — x + | is the second Bernoulli polynomial. 
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the theory deformed by a mass squared a, and then we make a into a dynamical field. The 
mass-deformed theory has the action (4.1), deformed by 



5S 



d 6 xa 



(4.35) 



The exact propagator (4.5) is now 
(W0t(_<?). 



~ 2 ,l y ^)H\p -q)= 1 5 Jt (2nf5\p - q) . (4.36) 

P 1 + (7 — 2^B P — ^B,a 

The computation of the bootstrap equation for the self-energy in section 4.1 generalizes to 
this case with minor modifications. The right-hand side of equation (4.18), given by the sum 
of 1PI diagrams, is given by the old one with the replacement T,b — > £b,<t- If we replace the 
definition (4.6) with S_b j(7 = —(3~ 2 ^, the right-hand side is unchanged. In the left-hand side 
of (4.18) we now have = £b, ct + a = —f3~ 2 (fj 1 2 3 — a), where a = j3 2 a. As a result, the 
final bootstrap equation (4.23) becomes 

\X\Hb 1 



±JfS B -a 



Li 2 (e-^- m|A| ) -Li 2 (e-^ s " 



■iri\X\\ 



(4.37) 



2 2iri 

The free energy written in terms of the self-energy is also slightly modified, and becomes 



,F B „ = W 2 r *, / *>, £ [.o g { f - E B ,„) + t^-^ 



2 Sr., + a 



(4.38) 



The extra contribution in the last term can be easily computed analogously to (4.32). The 
final result for the free energy is given by 



NV 2 T 2 f^3 2(/4-a) 3 / 2 , 1 



2?r 



3|A| 



7T*|Aj 



dyy 



II B 



Li 2 ( e -y+™W) - Li 2 ( e -^l A l 



(4.39) 



Let us now carry out the path integral over a to obtain the critical theory. The saddle 
point equation for a, derived from (4.39) (using (4.37)), is 



(4.40) 



and we find a = ii 2 B . Plugging this back in, we find the free energy of the critical bosonic 
theory, 



crit. 
B 



NV 2 T 2 /4, c | 1 

2tt I 3 TTlX 



dyy 



Li 2 (e 



— y+%i\\ 



Li 2 (e 



—y—niX 



where hb,c is given by (4.37) 



= A//R c + 



7T2 



(4.41) 



(4.42) 
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0.0 0.2 0.4 0.6 0.8 1.0 

X 

Figure 3: Thermal mass in the regular bosonic theory with Xq = (blue, lower curve) and 
critical theory (red, upper curve), with arbitrary normalization. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 4: Negative free energy in the regular bosonic theory with A6 = (blue, upper curve) 
and critical theory (red, lower curve), using arbitrary normalization. 
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These results are shown in figures 3 and 4. It is amusing to note that (4.42) is simply given 
by setting the \ib,c derivative of (4.41) to zero. 

One may check that in the limit A = these results agree with the thermal mass and free 
energy of the critical O(N) model computed in [37]. At A = 0, one finds from (4.42)-(4.41) 

(4.43) 

(4.44) 

5 7T 

which is indeed the result of [ 7] (up to a factor of 2 in the free energy, due to the fact that 
we work with U(N) rather than O(N)). 




5 Fermionic Theories 



In this section we consider, following [1], the theory of N Dirac fermions ip l , coupled to the 
U(N) Chern-Simons gauge field A^. The theory is defined by the action 

S = Sfermion + SCS> (5-1) 
Sfcrmion = J d 3 X frfD^ + U^lj) . (5.2) 

The regular theory is given by setting a = 0. In the critical theory (the Gross-Neveu model 
with Chern-Simons interactions) a plays the role of the primary scalar operator, and we must 
perform the path integral over it. 

The calculation is similar to the one done in [ ] . Our conventions for the gamma matrices 
are 7^ = a^, [i = 1, 2, 3, where are the Pauli matrices. In the thermal theory the fermions 
have anti-periodic boundary conditions on the thermal S 1 (before including the effects of the 
holonomy), so the momenta in this direction are quantized as ps = ^(n + ^), n £ Z. 



5.1 Exact propagator 

Denote the exact propagator as 



{2ir) 6 8 6 {p-q). 



(5.3) 



The notation is the same as for the bosonic theories, with Tip being the sum of 1PI diagrams. 
The bootstrap equation for Ejr is shown in figure 5. The 2x2 matrix Y^p can be uniquely 
expanded as the linear combination 



(5.4) 
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p 



a 



Figure 5: Bootstrap equation for Hp. A filled circle denotes the exact fermion propagator. 



Let us also define T,j = — a. The diagram in the bootstrap equation is given by 



10 



k jm ki J q (q-pY 



■7 



i 



7 



[3| 



iq^ + a-T F (q) / mk 
1 



7 



1+] 



27ri r m f, 1 

J (q-p) + iq^ + a - E F (q) 



1+] 



2vri r m r , 1 
— ^iTrjv / D 3 g 



1 



(9-p) + (g-S F )2 + S2 



7 



[3| 



i(g-E F ) M y* + EiJ 



7 



1+] 



(5.5) 



where in the second line we used (4.9). The right-hand side is proportional to the unit matrix 
in color space, so we may write Tp(p)ji = T,p(p)5ji. Using the relations ^PI^^W = —25^1 
and [7 3 ,7 + ] = 27 + , the bootstrap equation becomes 

v / \ 4 "t f-rf 1 £;h)7 + - i(g - Sf) + J 



~ 1/2 /" d 2 g 
-1/2 7 (2vr) 2 (g-p) + n 



I 1 ~ S7(g)7 + -i(g-^)+J 
^ {q-Tpf + t} 



We see that S_ = S3 = 0, and that Tip is independent of p 3 . Using also rotation invariance 
in the p 1 — p 2 plane, we may therefore write 



£j(p) = /(A, Pps, cr)p s , S+(p) = g(A, f3p s ,a)p + 
where a = /3<r. Let us now compute the sum in (5.6). As before, let x 



^ 00 



1 



1 



Z 3 n f^o (? " Sf) 2 + S 2 (g) /? V ?1 + (9 " Sf) 2 + S 2 

^2 

V [2vr (n+i-|A|n)] 2 + x2(l- 5 + /2) 



P 1 
2xVl- 9 (i) + f(x) 



Re 



tanh 



Vl- 5 (s) + / 2 (x) + 27ri|Ah 



(5.7) 



(5.8) 



10 Let 7 l Q U 7 l/31 = 7° ! J 47' 3 - 7 /3 A 7 a . 
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The bootstrap equation is now (comparing separately the coefficients of j + and /) 



£+(p) = 2vrA/3 [ du [ 

J-l/2 J 

E/0) = 2vrA/3 



d 2 q t!(q) 



-1/2 
1/2 



du 



1/2 



(2tt) 2 (g-p)- 

d 2 Q Q + 



(2tt) 2 



G(x) 
G(x) 



Let us differentiate with respect to p , using (4.14). We find 



d y g + g(y) = -yfG(y) , 



f(y) 



-G(y) 



Combining these, we find that 

d y (y 2 g-y 2 f 2 ) = 



g(X,y)-f 2 (\,y) 



/4(A) 



y 



(5.9) 
(5.10) 

(5.11) 
(5.12) 

(5.13) 



where [ip is an integration constant that may depend on the coupling. As we will see, \xf 
plays the role of the fermion's thermal mass. Returning to (5.10), we can write it as 



yf(y) + fa 



A f 1 ' 2 



2tt 
A 



du 



1/2 



dx xG(x) 



2tt 



dO n 



(q-p)- 



1/2 



du 

-1/2 

AA' - 2A 



dx 



Re 



X 2 + H 2 F 



tanh 



1/2 



du Re log 



-1/2 



2 cosh 




(5.14) 



Here again A' = (3 A (A being the cutoff). In the second line we computed the angular integral, 
which evaluates to 27r0(g s — p s ), and substituted G(x). In the third line we took A — > oo, 
keeping diverging terms. We subtract the linear divergence with a mass counter-term. For 
a trivial holonomy {u = 0) we find agreement with [1]. For our holonomy the remaining 
integral can be computed, and we find 

1 



yf{y) + Z = -\Jy 2 + y 2 F + 



Til 



Lio 



c.c. 



(5.15) 



To determine the integration constant uf in (5.13) (which we choose to satisfy \if > 0), we 
demand that g is regular in the limit y —> ( 
this limit, unless [If satisfies the equation 



demand that g is regular in the limit y — > 0. The expression f 2 — ^ has ay 2 divergence in 



±//f(A) = d + \\if + 



1 



1X1 



Li 2 



c.c. 



(5.16) 



It turns out that this equation always has a unique solution (with some choice of sign). 
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Let us summarize. The exact fermion propagator in the thermal theory is 



ty{p)d{-q)j) = r— . 1 t , T 5 jt (27T) 3 8 3 (p-q) 
iptfi* - igp+j+ - fp a I 

= — W+t*i4 ji( ^ (p ~ q) ' ( ^ 

where g is given by (5.13), / by (5.15), and \x F obeys (5.16). Notice that T\fi F \ is the thermal 
mass. 



5.2 Free energy 



For trivial holonomy, it was shown in [1] that the free energy F has a simple expression in 
terms of the self-energy. In our case the expression is similar, with the same modifications 
that were needed in the scalar case (4.25): 



-NVo 



1/2 
-1/2 



du 



d 2 



(2tt)* 



E ^ 



log(iq^-Z F (q))+lz F (q)- ' 



2 ' 11 iq^ - E F (q) J 
(5.18) 



where T, F = iS^7^ + S/J = — cr/. Let us first compute the trace over the fermionic 
indices : 



2 - Ej?(g) 



log 
log 



detfer.fig^ - S F (g) 



(g - £ F ) 2 + S| 



+ ^ Tr fer. 



Sf(9) 



-(q-X F f-£j t 



+ 



1 (5 " V' Ws ~ 2/(79. 



2 (9-£ F )2 + s2 
Notice that 

(g - S F ) 2 + E 2 = gf + {q- S F ) 2 + E 2 = q 2 + g 2 (l - g + f) = q 2 + 0~ 2 p F , 
where we used (5.13). The free energy can therefore be written as 

1/2 r fin J2. r 1 ( 5 - 2/ 2 )g 2 - 2faq s 



PF F 



-NV 2 



du 



1/2 



d z q 



E 



log (g 2 + /T >f) + 



2 q 2 + 13-^1 



(5.19) 



(5.20) 



(5.21) 



We first consider the log term. The sum over momentum in the compact direction diverges, 
and we regulate it as before by subtracting the contribution from a massive fermion with 
mass /3 _1 M: 



f d 2 q A 

7(2^ £ log 



2vT/3 2 



q 2 + P~ 2 M 2 
A' 

da; x log 



d 2 c 



(2tt): 



bg n 



1 + 



(27r(n + |-|A|u))^ 



2 cosh( Jx 2 + + 2 cos(27r|A|u) 



(mf -> Af) 



- (up -> Af) 
(5.22) 
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The remaining radial integral is the same as that of (4.26) if we replace fiB by /J>F and flip 
the sign of cos(27r|A|«). After subtracting the divergences as we did in the bosonic case, we 
find the result 



1 



2vr/? 2 



/4 



POO 

+ / dyylog(l + e- 2?/ + 2e" ?/ cos(27r|A|n)) 
Jfi F 



(5.23) 



1 /2 

We now compute the integral , 2 du of this expression. The result can be obtained from 
(4.30) by analytically continuing y — > y + ni, and we find 



27T/3 2 



t4 , i 

3 



dyy 



l-'F 



Li 2 (-e-^l A l)-Li 2 (-e- 



y-m\X\\ 



(5.24) 



Let us now consider the second term in (5.21). The sum was already computed in (5.8), and 
we find that 
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(5.25) 



Notice that from (5.15) we have 
d(xf(x)) x 



dx 



x 2 + /4 



log 
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X 2 + jl 2 p — TTlX 
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(5.26) 



Making the replacement v = xf(x) in (5.25) therefore allows us to solve the integral, with 
the result 
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(5.27) 



Combining (5.24) and (5.27) into (5.21), we find that the free energy is given by 
NV 2 U% 
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(5.28) 



The choices of signs here and in the previous equation depend on the sign for which there is 
a solution to (5.16). The thermal mass and free energy in the regular theory are obtained by 
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setting a = in (5.16) and in (5.28). We can always rewrite (5.28) as 



p F 2vr/3 2 \ 6A 6A TriX 



,2 



^Li 2 (_ e -^+- A ) + J dyyU 2 (- e ^+^) - c.c. 



(5.29) 



We can now perform a first test of the duality between the fermionic and bosonic theories. 
The duality implies that the expressions we found for a = should map to the analogous 
expressions (4.42) and (4.41) of the critical bosonic theory. One can verify that this is indeed 
the case, using the duality map 

|A 6 | + |A/| = 1, tt^tvt, sign(A b ) = -sign(A/) , (5.30) 
\*b\ |A/| 

where Nf,, Xb (Nf, A/) are the parameters of the critical bosonic (regular fermionic) theory. 
In particular, we see that in the A — > 1 limit the thermal mass fj,F of the regular fermion 
theory goes to the thermal mass (4.43) of the ordinary critical boson, and the fermion free 
energy goes to 

PF F ^ -NV 2 T\l - X)\^ ~ -k YM V 2 T^ C -^ , (5.31) 

5 7T 5 7T 

which is the free energy of the bosonic critical C/(/cym) model [37]. 



5.3 The critical theory 

Consider next the critical fermionic theory. Like the regular bosonic theory, this has (at large 
N) an exactly marginal triple-trace deformation, that we should add also on the fermionic 
side in order to map to the bosonic theory with general values of A6- At large N, this theory 
can be defined by deforming the regular theory by <5£f crm i on = cj^^ + ^AgC" 3 , computing the 
effective action for a, and finally computing the path integral over a (which at large N is 
equivalent to extremizing the effective action with respect to a). The free energy (5.28) is the 
effective action for a in the thermal theory, and turning on the Ag deformation corresponds 
to adding a term 

6(l3F F ) = ^4d*. (5.32) 
Extremizing (3Fp with respect to a, we obtain 

a = -sign(A)/i F 5 , 9 = , , (5-33) 



1 - 2vrAA 



where the sign of a was chosen such that there is a solution of the thermal mass equation 
with up > for any A and Ag such that 27rAAg < 1. The dominant solution (with the lowest 
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free energy) has the sign in (5.16) equal to the sign of A. The free energy in the critical 
fermionic theory is then given by (with an '/' subscript on fermionic couplings) 
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(5.34) 



where fj,p c obeys 



(1 - |A/| + \gf\)nF,c 



1 

TT'l 



Li2 



C.C. 



(5.35) 



For A = A6 = the free energy of the critical fermion theory is the same (at large N) as that 
of the free fermion theory (see [38, 39] and references therein), but this is not true for finite 
couplings. 

Let us now compare this result with the free energy of the regular bosonic theory. Up to 
normalization, a is dual to the scalar operator Jo = 4> <t> in the regular bosonic theory, and 
the a 3 interaction is dual to the bosonic triple-trace deformation <5£boson = ^^(Jo) 3 - By 
matching the 3-point functions (aaa) in the critical-fermionic theory with (JoJoJo) in the 
regular bosonic theory, the mapping of the triple-trace couplings in the two theories is [( 
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(5.36) 



Under the duality map (5.30) and (5.36), we have that = Xl + — 



b + 8^ maps as 



\h\ 



2 

\9f\ 



(5.37) 



One can now easily check that under this map the free energy (5.34) and thermal mass (5.35) 
of the critical fermionic theory map to those of the regular bosonic theory, given by (4.33) 
and (4.23). This is true for any value of A and Xq on both sides, with the constraint 



xU 1 



Xl > 



-8TT 2 X 2 b 



(5.38) 



6 Theories with Bosons and Fermions 



In this section we generalize the computation of the free energy to massless theories with both 
fermions and scalars, with generic couplings. As a special case we will consider the theory 
with J\f = 2 supersymmetry (for a review, see [10]). As a test of the computation, we will 
verify the Seiberg-like duality in this theory that was found in [24] (generalizing the duality 
for non-chiral theories found in [23]). 
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5Y*f = ^ > 5Y*b = 

(a) (b) 

Figure 6: Additional diagrams due to the A4 vertex, in the bootstrap equations for the exact 
propagators. 

The theory we consider includes a single complex scalar and a single Dirac fermion, both 
in the fundamental representation of U(N), coupled to a gauge field with Chern-Simons 
interactions. It is defined by the action 

S = S C S + Scalar + Sfermion + j d 3 X ^ (lp1p)((j)^ (j)) , (6.1) 

where the first three terms are defined in (3.1), (4.2), and (5.2), and we take the fermion and 
scalar to be massless for simplicity (a = 0). 11 

6.1 Exact propagators 

The scalar and fermion propagators are given in (4.5) and (5.3), respectively. There is one 
additional diagram for each bootstrap equation, shown in figure 6. These diagrams are given 
by 

•V3 r rfir 



(a) = - w / 1/2 du Iwf \ ? F^U) 



(6.2) 



(b) = M * IZ dU lw?~^ • (6 - 3) 

Using the same arguments as in the separate bosonic and fermionic theories, we find again 
that £_ = S3 = 0, and that 

Sb(p) = -/T 2 m!(A) , (6.4) 

X I (p) = f(y,X)Ps, (6.5) 

s +(p) = g{y,^)p+, (6.6) 

where 

9-f = -^. (6.7) 

yZ 

11 One can write down additional marginal interactions of the form and ((?/>0)(V></>) + ex.), but 
they will not affect the free energy in the large N limit. 
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The diagrams can be computed using the same techniques we used above. Let us define 

1 r 



h B = \X\fi B + 



IT! 
1 



Li 2 (e-^-^ |A| ) -c. 



c. 



hp = \X\fiF H : Li 2 (— e 



TTl 



-fl F -TTl\\\\ 



The integrals in (6.2) and (6.3) then evaluate to 
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(6.11) 



To compute the last integral in (6.11) we first determined f(y), which using the results of 
the previous section and (6.10) evaluates to 



yf(y) = -XJy 2 + fi 2 F + 



sign(A) 



Lis 



-y/y 2 +ti 2 F +7Ti\\\ 



c.c. 
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(6.12) 



Adding the contributions (6.2) and (6.3) to the bootstrap equations (4.18) and (5.15), we 
find the corrected equations for the thermal masses, 

Ml = (±h B ) + ±h F (h F - ^h B ) , (6.13) 



Hf = hp 



X^_ 
4vrA 



h B . 



(6.14) 



In the second line we chose the branch that gives a solution with positive fj, F . 
6.2 Free energy 

The free energy receives separate contributions from the scalars and fermions, as well as 
contributions from the A4 interaction. It is given by [4] 



F = F B + F F + F BF 



(6.15) 
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where F B is given (in terms of the self-energy) in (4.25), Fp in (5.18), and 
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(6.16) 



In the second line we used the integrals (6.10) and (6.11). 

The bosonic contribution Fb to the free energy is given by (4.34), and the fermionic 
contribution can be computed as in section 5.2. The only difference is in equation (5.27), 
where we must account for the corrected thermal mass equation when computing the integral. 
We find that 



1/2 



du 



d 2 



1/2 7(2tt)2/3 



. 00 
R 



(g~2f 2 )q 2 



1 



P-*H% 8tt|A|/33 



Ai 



4vrA 



h B 



+ 



1 



A. 

4ttA 



h B 



(6.17) 



The resulting fermionic contribution to the free energy is 

V 2 N j \X\t4 m| [i\h F 1 
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(6.18) 



Combining all the contributions to the free energy and using (6.8), (6.9) and (6.14) to 
simplify the expression, we obtain 
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(6.19) 



6.3 The supersymmetric case 

The theory with M = 2 supersymmetry is conjectured to be self-dual [ ], with the transfor- 
mation at large given by N — > \k\ — N and k —> —k. In this section we verify that the free 
energy is invariant under this duality. 

In this theory, the couplings are given in terms of A by 12 [ ] 

A 4 = 4vrA , A 6 = 24vr 2 A 2 , (6.20) 

12 The supersymmetric theory also contains other couplings which do not contribute at large N, as mentioned 
in footnote 11. 
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Figure 7: Free energy (times —A) in the N = 2 theory, in arbitrary normalization. 



so that A = 2|A|. The equations for the thermal masses become 

fi B = (h B - h F ) 2 , (i F = h F - h B ■ 

The thermal masses are therefore equal, {j, = fj, B = fip, with /i satisfying 

1 
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The free energy is given by 
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(6.21) 



(6.22) 



(6.23) 



Both the thermal mass and the free energy are invariant under the duality, which in terms 
of N and |A| takes |A| — > 1 — |A| with N/\X\ fixed. The free energy is plotted in figure 
7. It is straightforward to generalize our results also to theories with different amounts of 
supersymmetry, including the theories analyzed in [41, 4]. 



7 Summary of Results 

In this section we collect the main results obtained in this paper. 

We have computed the thermal free energy on R 2 for various large N vector models coupled 
to Chern-Simons gauge fields, working exactly in the 't Hooft coupling A = ^. In particular, 
we have explicitly tested the conjectured non-supersymmetric dualities ("3d bosonization" ) 
relating the regular /critical fermion coupled to Chern-Simons to the critical/regular scalar 
coupled to Chern-Simons [1, 8, 5]. Our calculations closely followed [1, 4, 20], with the 
difference that we included the effect of the holonomy around the thermal circle, as explained 
in Section 2. 
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Figure 8: Plot of (—A) times the free energy, (— \Fp), (blue solid line) and A times the stress- 
tensor 2-point function X(TT) ~ sin ^ x (red dashed line) for the regular fermion theory. They 
look very similar but they have a different functional form. 



In all cases we find that the free energy goes to zero at A = 1, in accordance with the 
behavior of the exact (TT) correlator [5, 6]. In fact, since the thermal free energy F and 
(TT) are both a measure of the number of degrees of freedom of the theory, it is interesting 
to compare them. Amusingly, we find that their behavior as a function of A is rather similar. 
As an example, figure 8 shows a plot of (-XF) and \{TT) for the regular fermion theory, 
normalized to agree at weak coupling. Note that it is possible to smoothly continue our 
results beyond A = 1. In fact, the functions XF and X(TT) in the regular fermion theory 
(and the critical bosonic theory) are both periodic under A — > X + 2 and odd under A —¥ —A. 

Our results are most simply expressed in terms of the following functions 
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(7.1) 

(7.2) 
(7.3) 



7.1 Regular fermion and critical boson results 

For the regular fermion theory (a massless fermion coupled to CS), the exact planar free 
energy takes the form 

NV 2 T 3 



X 



-F{Hf, A - sign(A)) 



where the thermal mass fip is determined by the equation 

[i F = h F (fJ.F,^) ■ 



(7.4) 
(7.5) 
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Note that here (unlike our discussion above) we allow hf also to be negative, and the equations 
are still valid in this case. In fact, we have hp{—pLp, —A) = —hp(fiF, A), and the free energy 
is even under A — >■ —A. For the critical boson coupled to CS, we found 

Fr- = ^^H»B,cA), (7.6) 
where ^b,c is given by (again allowing also negative values) 

h B (fiB,c,*) = 0, (7.7) 

It is now obvious that if the thermal masses map to each other under the duality, so will 
the free energies. Note that y = Jcym + ^ r sign(/cYM) changes sign under the duality It is 
easy to see that the thermal masses are indeed related by |//f(A)| = |//b jC (A — sign(A))|, in 
agreement with the expected duality. 13 The inclusion of the non-trivial holonomy was crucial 
in order to obtain this. 

Our theories have a global U(l) symmetry, which acts on the scalars/fermions in a natural 
way (the equations of motion of the gauge field identify this symmetry with the topological 
symmetry of the three dimensional £7(1) gauge field). It is simple to generalize our analysis 
by adding a chemical potential for this gauge field, as done (without taking the holonomy 
into account) for the regular fermion theory in [42]. In the bosonic theories we find that the 
effect of adding a chemical potential fijj^ is simply to replace every Li2(e x ) appearing in our 
analysis by ^(Li2(e x+>lu W) + hi2(e x ~ >lu W)). The chemical potential maps to itself under the 
duality. 

For comparison to results that previously appeared in the literature, one may also look at 
the perturbative expansion of these results around A = 0. The first few orders of the thermal 
mass and free energy for the regular fermion theory read 

H F = 2 log(2) A + (log 2 (2) - ^) A 3 + 0(A 5 ) , 



F F = -NV 2 T S 



3C(3) _ 7r 2 log(2) + 41og 3 (2) 2 ^ 
4tt 6vr { J 



(7.1 



One may verify that these differ from the results in [1] already at order A 2 . Analogously, one 
may derive the perturbative expansion of the critical boson free energy. Since the thermal 
mass ^b jC (A = 0) = 2 log f v ^ +1 ^ is non-zero, the perturbative expansion is analytic. By 
the duality map, this also means that the expansion of the fermion free energy (and thermal 
mass) close to A = 1 is analytic, unlike the result found in [I]. 



Note that the function J-(n, A) is even in fj,, so the thermal free energy is insensitive to the sign of fj,. 
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7.2 Regular boson and critical fermion results 



In this case, the situation is slightly more complicated due to the possibility of adding a 
sextic coupling ^^(Jq) 3 , where Jo is the scalar/fermion bilinear. In the large N limit, this 
coupling is exactly marginal [2] and we have computed the free energies keeping A6 as a free 
parameter. The result for the regular boson turns out to be 

NV?T 3 

F B = —±—F{^\), (7.9) 

where the thermal mass is determined by 

2A 

— w-hb = h,B(fJt>B, A) , (7.10) 
A 

with A 2 = A 2 + A. 

For the critical fermion theory, we found that the free energy is given by 

F crit. = _^l HflFcj A _ sign(A)) ( 

and the thermal mass is given by 

(l + g)VF,c = h F (vF,cA), (7.12) 



where g = 1/ y 1 — 27rAAg . 

In this case, the duality requires a non-trivial map between the \q couplings, given in 
equation (5.36), which can be derived independently by matching the scalar three-point 
functions in the two theories [6]. For example, the critical fermionic theory with Ag = 
(g = 1) maps under A •<->■ A — sign(A) to the regular bosonic theory with A6 = 24-7r 2 A 2 
(A = 2|A|), as can be easily seen from the expressions above. 

7.3 Theories with bosons and fermions 

Let us now summarize the results for the theory with one scalar and one fermion, discussed 
in section 6. This theory contains a boson- fermion coupling ^-(^ip)^^) in addition to the 
scalar sextic coupling g^(0^) 3 . We found that the free energy in this theory is 

NV?T 3 

F = (H^b, A) - T&F, A - sign(A))) , (7.13) 



where 
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HF = h F --^h B , (7.15) 
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and A is denned as in the regular boson theory. 

The M = 2 supersymmetric Chern-Simons theory with one chiral superfield is a particular 
case of the above theory, obtained by setting A = 2|A| and A4 = 4tt\. In this case we obtain 

\i = li f = [i B = hp - h B ■ (7-16) 

One finds that ll is self-dual under the level-rank duality transformation, consistent with the 
expected Seiberg-like duality of this theory [24]. 

One would expect that the theory with one scalar and one fermion is self-dual (at infinite 
N) also away from the N = 2 supersymmetric point, for general values of A4 and Xq, since 
we can deform the two theories away from the supersymmetric point, and the deformation is 
exactly marginal for infinite N (but not for finite N). Let us now describe how this duality 
works in the large N limit. 

The single-trace operators of this theory consist of the operators of the regular boson 
theory, which we will denote by J% (s = 0, 1, • • • ), and those of the regular fermion theory, 
denoted by j/. In addition, there is a tower of half-integer spin operators obtained by 
contracting the scalar with the fermion. 

It is easy to generalize the computations of planar correlators of integer spin operators, 
performed in [ , ], to this theory. One finds that the results are consistent with conformal 
invariance for any A4 and Xq. In particular, by computing the exact planar 2-point functions of 
the operators (ftfyfyip) and ($(j)) 3 we find that these operators have dimension 3 + 0(l/iV), 
consistent with the expectation that these deformations are exactly marginal in the planar 
limit 14 . 

The planar 3-point functions of the operators J% or j/ with s > 1 are independent of A4 
and A6, giving the same results as in the regular boson or fermion theories. Therefore the 
correlators of the currents Jg map to those of j/ under the level-rank transformation. The 
A4 deformation affects planar 3-point functions with insertions of the scalar operators Jq and 
Jq, while X$ affects only the 3-point function of Jq- 15 By computing 3-point functions with 
T b 



insertions of Jq and j/ we find the duality map of A4 and Xq to be 



-4 - ^ - ^ , (7.17) 



47rA X4 
4A 2 x\ 



l + xi-—\. (7.18) 



Note that the above transformations are consistent with the self-duality of the N = 2 theory 

A 2 
4A 2 



in which X4 = ^ = 1. 



14 This fact is also consistent with the perturbative computations of the beta- functions, done in [43]. 
15 Note that Jq and Jq have different dimensions in this theory (1 + 0(1/N) and 2 + 0(1/N) respectively), 
so they must be self-dual. In particular we find that Jq — > — AA4 Jq, up to some positive multiplicative factor. 
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One can now verify that the free energy (7.13) maps to itself under the above transfor- 
mations for any values of A4 and A6, consistent with the expected self-duality of this theory. 
Note that to complete the duality map one has to determine also the mapping of the half- 
integer spin operators, and include the effects of the additional double-trace deformations of 
this theory, constructed from the spin 1/2 operators (i/>0) and (<$ip) (see footnote 11). While 
the free energy does not depend on those deformations, correlators of the spin 1/2 operators 
will. Moreover, one has to prove that all the above deformations are exactly marginal more 
carefully. We leave a more detailed analysis of these issues to future work. 



7.4 Massive theories 



While in this paper we have concentrated on the conformal theories, we should note that our 
results also allow us to readily obtain the thermal free energy for massive fermions or scalars 
coupled to a Chern-Simons gauge field. Studying these theories in detail is beyond the scope 
of this paper, but let us make a few brief comments about them. 

For example, the regular fermion free energy derived in (5.28), with [ip obeying (5.16), 
describes a massive fermion theory, where the mass parameter appearing in the Lagrangian 
is mp = a/ f3. For convenience we rewrite here the two equations defining fhf = (3m f, 
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(7.20) 



As above, we allow also negative values of hf] there is always a solution for either positive 
or negative [ip. 

This theory should be related by the duality to the critical boson theory, with a mass 
deformation mj, proportional to mj. In this theory we continue to get equation (4.37), but 
now we need to add an extra term to the free energy in (4.39) of the form — inside the 
curly brackets of (4.39) (where again rhb = (3rrib). The equation for a is now modified into 

4 



Substituting this value of a in (4.37) and (4.39), we obtain (using here /ib > 0) 



-Xrhb = —X/j-b 
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(7.22) 
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(7.23) 



In writing (7.23) we added to the Lagrangian density a constant term proportional to NX 2 mf 
(this modifies only the second term, and does not affect the physics), in order to match with 
the fermionic theory, using the identifications 

\Hb\ = \hf\, \b = \f — sign(A/), X b m b = -m f . (7.24) 

The last relation between the mass parameters follows (up to a sign) from properly normal- 
izing the associated spin zero operators, using their two-point functions computed in [5, 6]. 

Actually, this is not precise, since equation (7.22) does not always have a solution (recall 
> 0). When it does not, it turns out that there is (at least at zero temperature [5]) a 
different configuration in which one of the scalars obtains an expectation value, breaking the 
U(N) gauge symmetry to U(N — 1). The duality implies that analyzing this configuration 
(using the expectation value that minimizes the free energy) should lead to exactly the same 
equations that we had above, just for negative values of (j,b, in order to be consistent with 
the matching with the fermionic theory for these values of the mass. It would be interesting 
to confirm this. 

It is interesting to take the zero temperature limit of the above relations. In the zero 
temperature limit rh and \i become large, so the dilogarithms go to zero if the argument 
of their exponent is negative (otherwise we can map to this situation, using the identity 
Li2(x) +Li2(l/x) = — ^ — ^(log(-x)) 2 ). For the two theories described above, we find the 
simple results (assuming for simplicity 0<A&<1, — 1 < A/ < 0) 

T = ~T =mb = ~ (l+ f \ f ) ' m b >0,m f <0, (7.25) 

HB A 6 m 6 m f 

y = -y = (2^ = -(T3A7) ' ™ b <0,m f >0. (7.26) 

Note that /if,b//3 is the position of the pole of the dressed fermion/scalar propagator in 
these theories. This is not a gauge-invariant quantity. However, at least for weak coupling, 
we expect that the position of this pole is close to the actual physical mass of the lightest 
(charged) excitation. In the massive theory, the low energy effective description is by a 
topological pure Chern-Simons theory. We have massive excitations whose worldlines act 
as Wilson line sources for this Chern-Simons theory, with vanishing forces between them at 
leading order in the large N limit. 

In the fermionic theory for small A/ the excitation mass is of order nif. However, it seems 
strange that when A/ is not small, we have a significant difference between the two choices of 
signs for nif] the low-energy theories in the two cases differ just by shifting k by one (which is 
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negligible in the large N limit). However, we can understand the origin of this difference by 
thinking about the critical bosonic theory. In this theory we see that for small A;, the physical 
mass of the excitations is the same as for mj, > 0, but is a much smaller scale of order 
AfeTOfc for negative m&. As mentioned in [5], this is because for one sign the theory is in a phase 
with unbroken U(N), while for the other sign the U(N) symmetry is spontaneously broken 
to U(N — 1) by a scalar vacuum expectation value. Since the U(N) symmetry is gauged, 
instead of having Nambu-Goldstone bosons, we have light massive gauge bosons with a mass 
of order XbiJib (for At, <C 1). 

There is also another way to obtain massive theories. Instead of a massive deformation, 
in some cases we have (for zero temperature in the regular bosonic and critical fermionic 
theories) a moduli space along which we can spontaneously break the conformal symmetry 
by a vacuum expectation value for the gauge-invariant scalar operator. Without the Chern- 
Simons coupling this possibility was analyzed in [44, 45]. Usually the only solution to our 
equations at zero temperature is fj, = 0. However, in our regular bosonic theories, we see 
that when A = 2 (namely A 2 + = 4), there is a solution to the "gap equation" (7.10) for 
the thermal mass for every negative value of ( m the zero temperature limit). This is the 
manifestation of the opening of a new moduli space of vacua for this specific value of A in 
our analysis. Along this moduli space the conformal symmetry is spontaneously broken, and 
there is a massless Nambu-Goldstone boson. Of course, for finite temperature this moduli 
space is lifted, and we also do not expect it to be present for finite values of N. 
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